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The POWHEG method applied to top pair production and 

decays at the ILC 

Oluseyi Latunde-Dada * 

Cavendish Laboratory, University of Cambridge, 
JJ Thomson Avenue, Cambridge CB3 OHE, U.K. 

We study the effects of gluon radiation in top pair production and their decays for e + e~ 
annihilation at the ILC. To achieve this we apply the POWHEG method and interface 
our results to the Monte Carlo event generator Herwig++. We consider a center-of- 
mass energy of yfs = 500 GeV and compare decay correlations and bottom quark 
distributions before hadronization. 
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1 Introduction 

In Tabic [TJ we have highlighted some differences between matrix element (ME) and partem 
shower (PS) generators and have labelled as (M) or (D) those attributes we consider merits or 
drawbacks respectively. Note that most PS generators attempt to include NLO corrections 



PS generators 


ME generators 


Resums leading logarithmic contributions 
to all orders (M) 


Can only go up to N~ B LO (D) 


High multiplicity hadrons 
in the final state (M) 


Low multiplicity partonic 
final states (D) 


Works well in regions 

of low relative pt (M & D) 


Works well in regions 
of high relative p T (M & D) 


Total rate is accurate to LO (D) 


Total rate is accurate to N( >0 )LO (M) 



Table 1: Differences between PS and ME generators 

via a method called the matrix element correction which corrects the hardest shower emission 
so far to the exact matrix element and populates the high px regions according to the NLO 
cross-section. However, the total rate is still only accurate to LO and virtual corrections are 
not fully taken care of. 

1.1 Getting the best of both worlds at NLO 

The Positive Weighted Hardest Emission Generation (POWHEG) method |H [3] aims to solve 
this problem. It 

1. generates total rates accurate to NLO, 
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2. treats hard emissions as in ME generators, 



3. treats soft and collinear emissions as in PS generators, 

4. and generates a set of fully exclusive events which can be interfaced with a hadroniza- 
tion model. 

The POWHEG method achieves this by generating the hardest emission in the shower first 
to NLO accuracy using a modified Sudakov form factor. For angular ordered showers like 
Herwig++, it also includes a truncated shower of soft, wide angled emissions from the hard 
scale to the scale of the hardest emission. This maintains the correct soft emission pattern. 
This is illustrated in Figure [TTT1 It then showers the resulting partons subject to a px veto 



Hard process 





Hardest emission 



Truncated shower 



Figure 1: POWHEG emissions 

to ensure that no harder emissions are generated. Unlike MC@NL0 [1] , it is independent of the 
PS generator used and all events have positive weight. In this talk [I], we will focus on the 
description and applications of the method in conjunction with the PS generator, Herwig++ 
0. 

1.2 The partem shower hardest emission cross-section 

For a single parton, the cross-section for the hardest emission with transverse momentum 
Pt is given by 

Av(0) + A V (PT)S- Pd *% 



da = dav 



2tt 



(1) 



where P is the splitting function for the hardest emission and Av(pr) is the Sudakov form 
factor for no emissions with fcx > pt which is given by 



Av(pt) = exp 



(2) 



The cross-section (JTJ expanded to order as gives 



da = dav 
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(3) 



The POWHEG method aims to substitute © with the exact NLO result within the parton 
shower. 

1.3 Correcting to the exact NLO cross-section 

The exact NLO cross-section can be written as 



rfoNLO = ^ ct b + day + c?ctr = daB + day + da^drM 



(4) 
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Adding and subtracting da B J s dr(M — C) we get 

dam,o = da B + day + da B / dr(M — C) + da B drM — da B / dr(M — C) . (5) 

Js Js 

where C is a counter-term and 5 is the subtraction region. This can be rearranged to give 



do nlo = day + d<r B J dr(M - C) + da B jl — J drA4^ + Mdr 



(6) 



with day/ = day + da B L drC now finite. Comparing ^ with ^ above, we can write down 
an analog of |T]) as 

da NLO = da B [A NLO (0) + A NLO (p T )Mdr] (7) 



where 



da B = da B + day + da B / dr(M. — C) 

Js 

J MO(k T -p T ) 



Anlo(pt) = cxp 



(8) 



Note that in defining da B , we have neglected terms of higher order than as and if it is 
negative, perturbation theory has broken down. 

1.4 POWHEG formalism and applications 

With 

dcr NLO = da B [A NLO (0) + A NL oOr).McZr-] , (9) 
the PDWHEG method can be applied by 

1. generating the pt of the hardest emission and its emission variables r, according to 
the term in square brackets using well known Monte Carlo techniques, 

2. distributing the underlying Born variables according to da B , (This defines the event 
weight and since it is always positive definite, all event weights are positive) 

3. for angular ordered showers, implementing a truncated shower of soft emissions be- 
tween the hard scale and the scale of the hardest emission, 

4. and finally showering the resulting partons as in a PS generator subject to a veto. 

The POWHEG method has been applied successfully to the following processes Z pair 
hadroproduction [B], heavy flavour production [7|, e + e~ annihilation to hadrons [BJ, Drell- 
Yan vector boson production [9l [10] and Higgs boson production via gluon fusion [11] . 

2 Top-pair production and decay at the ILC 

The application to top pair production and decay at the ILC takes the following into account. 

1. Spin correlations are taken into account in the matrix elements, M. for the production 
and decays of the top pairs. 
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Figure 2: Top-pair production and decay 



2. The narrow width approximation is applied so that production and decay interference 
can be neglected. This independence enables us to apply the method in separate 
frames: the lab frame for production and the top rest frame for its decay. 

3. In the lab frame, the transverse momentum fc^ is defined relative to the original t — t 
axis whilst in the top rest frame it is relative to the original b — W axis. 

4. The scale range available for production emissions (as log(-y/s/mt)) and is much less 
than the range available for decay emissions (sa log(m t /mb)). 

5. There are two different sources of the decay emissions: one from the top quark before 
it decays and the other from the b quark after the decay. Hence, there arc three 
different regions for truncated emissions labelled g(tr) in Figure [21 These are before 
the hardest emission in the production, from the top quark before it decays and before 
the hardest emission from the b quark. 

Further discussion of the method and its application can be found in |12j . Setting y/s = 500 
GeV and m t = 175 GeV, we considered the following four cases with PQWHEG interfaced with 
Herwig++: 

1. Leading order (LO) with no PDWHEG emissions, 

2. Only PDWHEG emissions in the production process (Prod), 

3. Only PDWHEG emissions in the decays of the tops (Dec), 

4. Both production and decay emissions allowed (Prod + Dec). 

For the two different e + e~ initial polarizations, we investigated the correlations between the 
decay products (we consider leptonic decays only for the W bosons) and the momentum 
distributions of the b quarks before hadronization. A selection of the plots obtained are 
presented in Figure 3. 
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Angle between the charged leptons from the top and anti— top decays 



Angle between the charged leptons from the top and anti— top decays 





Energy of the b quark (GeV) 



Energy of the b quark (GeV) 




Transverse momentum of b quark (GeV) 



Transverse momentum of b quark (GeV) 




Longitudinal momentum of b quark (GeV) 



Longitudinal momentum of b quark (GeV) 




Figure 3: Correlations and b quark momentum distributions. 
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3 Summary 



NLO improvements of parton showers are essential for near accurate predictions of angular 
correlations and momentum distributions at future colliders. The POWHEG method achieves 
this by distributing the hardest emission according to the NLO matrix element and yields 
events with positive weights. For angular ordered showers, the addition of a truncated shower 
is required. The method, though not very straightforward to apply, has demonstrated success 
in comparison with existing collider data. 

In this talk, we have extended this to top-pair production and decay at the ILC and 
made predictions for some distributions in comparison to leading order predictions. 
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